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Abstract 


In the hot and dense medium of a supernova (SN) core, the nucleon spins fluc¬ 
tuate so fast that the axial-vector neutrino opacity and the axion emissivity 
are expected to be significantly modified. Axions with ma ^ 10“^ eV are not 
excluded by SN 1987A. A substantial transfer of energy in neutrino-nucleon 
{nN) collisions is enabled which may alter the spectra of SN neutrinos relative 
to calculations where energy-conserving uN collisions had been assumed near 
the neutrinosphere. 
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Introduction. Nucleons interact with each other by a spin-dependent force which causes 
their spins to flip in a typical nucleon-nucleon (NN) collision. In the hot nuclear medium of 
a supernova (SN) core where p = gcm“^ and T = 10-60 MeV the spin-fluctuation 

rate To- is estimated to be at least of order T. Fluctuating spins act as a source for radiation 
such as neutrino pairs or axions which couple to the nucleon spin. Also, particles which 
scatter on nucleons by a spin-dependent force such as neutrinos or other nucleons can gain 
or lose energy if the target spin fluctuates. Further, fast spin fluctuations cause an averaging 
of the local spin density which reduces the effective vN or NN scattering cross section. To 
lowest order the physical essence of these processes is described by the Feynman graphs of 
Fig. §[] In a series of papers we have called attention to some of these phenomena 
and their consequences for the SN emission of neutrinos and axions. In the present note we 
provide a consistent account of this research and summarize its implications. 

The Spin-Density Structure Function. We focus on an isotropic, nonrelativistic, nonde¬ 
generate medium of baryon density hb, temperature T, and a single species of nucleons. The 
axial-vector neutrino interaction with this medium is governed by the dynamical spin-density 
structure function (e.g. Refs. 0,|]) 

4 r+ca 

S„{u},k) = - — / dte*‘^‘(cr(f,k) ■ cr(0,-k)) (1) 

oUb J-oo 

where cr[t, k) is the spatial Fourier transform at time t of the nucleon spin-density operator 
(t(x) = {x)T'ijj{x). Here 'ipi^x) is the nucleon held (a Pauli two-spinor) and r is a vector 
of Pauli matrices. The expectation value (...) is taken over a thermal ensemble so that 
detailed balance k) = Sa-{—uj, — k) is satished. 

The axial-current transition rate for a neutrino of four momentum (£i,ki) to { 62 , ^ 2 ) 
is then found to be \C\GYnB{‘i — cos 6 *) S'o-(a;, k) with 6 the scattering angle, (a;,k) = 
(£1 — £ 2 ,ki — k 2 ), and G-p the Fermi constant [^. The axial neutral-current coupling Ga 
to free p’s and n’s should differ somewhat from the charged-current value 1.26 because of 
the polarization of strange quarks in nucleons. Further, a suppression by ~ 20% of all axial 
couplings is likely in a bulk nuclear medium . 

We further neglect the momentum transfer from neutrinos to nonrelativistic nucleons 
(long-wavelength approximation) so that only Sa{oj) = limk^o S^ioj^ k) is used [^. After an 
angular integration the axial-current scattering cross section is 

daA _ ‘iG\Gl el 5'^(gi - 62 ) 
de 2 4:71 271 

This result is based on the interaction Hamiltonian (C'^Gf/2a/ 2) '0Ar7/i75'0iv " 01 / 7 /^( 1 - 75 ) V’i.. 
Axions a couple to nucleons by (CAr/2/a)where Gn is a model-dependent 
factor and fa the Peccei-Quinn scale P|. The axionic energy-loss rate per nucleon is 

~ = 7i%w/ dujuj^Sa{-uj). (3) 

ns (47r)Va do 


Wig.lJ) does not provide a complete lowest-order description of “inelastic nN scattering” because 
one needs to include wavefunction and vertex renormalization terms |18]. However, our use of a 
sum rule obviates the need to include these terms explicitly ||Tt 
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Within our approximations, Qa and a a depend on the same function Sa{uj). For a mixture 
of n and p the nuclear response is characterized by a linear combination of a small number 
of functions to allow for the isospin dependence of different processes. 

In a noninteracting medium the nucleon spins do not evolve, and so in Eq. (|l|) cr{t) = 
cr(0). Then the time integration yields = 27r5(a;). Indeed, if one ignores recoil effects 

the neutrino cannot gain or lose energy in a uN collision. Also, Qa = 0 in agreement with 
the fact that free nucleons do not emit radiation. 

Nucleon-nucleon collisions with a spin-dependent force cause a nontrivial evolution cr[t). 
Still, in doj S^{(jj) the factor gives us 27i6{t). Therefore, one hnds that the “sum” 
duj = (d/Sns) 27r(cr(0) •cr(0)) is independent of the time evolution of cr(f). There¬ 

fore, unless the NN interaction establishes spin-spin correlations which are absent in a non¬ 
interacting medium, one hnds the sum rule du S^(u) = 27r. 

In Ref. 0 the static spin-density structure function S'o-(k) = Sa-(ciJ,k) dcj/27r, and 
notably its long-wavelength limit k —0, has been used to study correlation effects for SN 
neutrino opacities. We believe that it is crucial to keep the dynamical structure because the 
cu-dependence is not essentially d(aj). Still, the results of P indicate a signihcant reduction 
of S'o-(k) due to NN interactions so that the partial pairing of spins may reduce du Sa{uj) 
below its free value 27r. 

Ansatz for S'o-(a;). Motivated by a classical bremsstrahlung calculation we assume a 
functional form S^{uj) = V„uj~‘^s{uj/T) h{ijj/T) where To- is the spin-huctuation rate, s{x) 
a dimensionless even function normalized such that s(0) = 1, and b{x) = 1 for x > 0 and 

e® for a: < 0 to satisfy detailed balance. If one models the NN interaction by a one-pion 

exchange (OPE) potential one hnds for a single species of nucleons 1^,0 

ho-,OPE = d-^/vr al -= 8.6 MeV pi 3 (d) 

where a.„ = (/2mAr/m^)^/d7r 15 with / 1.0 is the pion hne structure constant, pis = 

p/10^3gcm“^, Tio = T/lOMeV, and the vacuum value tun = OdOMeV has been used. 
In the following we shall use the classical result s{x) = 1. The detailed large-x behavior 
depends on the assumed short-range behavior of the NN interaction potential, but the OPE 
potential is too singular at short distances to yield a useful answer. 

For a; < Fa- the behavior of S„{uj) is dominated by multiple-scattering ehects. A proper 
calculation in this regime does not exist, but two of us (RS) have previously argued that a 
plausible representation is |l|J^ 

(5) 

^2 _j_ p^/4 

with F Fo-. Taking s(x) = 1 the sum rule is satished if one chooses F = Fa for Fa -C T, 
and F = Fa/2 for Fa ^ T. Assuming this form for Sa{uj) allows us to compute the variation 
of the thermally averaged axial-current scattering cross section (ua), and of Qa/’^s, with 
Fa (Fig. In a dilute medium (Fa -C T), {a a) is independent of density. The axion 
emission rate per nucleon Qa/ns is proportional to ub and thus to Fa. Such conditions may 
be relevant near the “neutrinosphere” where pis 0.1. However, for the conditions in a SN 
core Eq. (|D yields Fa,oPE/^ = 20-50 0. While we will argue that this perturbative result 
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likely is an overestimate of the true Fq- it still indicates that multiple-scattering effects may 
be expected to suppress both {a a) and Qa/ns- 

“Calibration” of Fq- at High Density. In the ansatz Eq. (j^) for it was implicitly 
assumed that Fg. increases linearly with density, as for Fq-ope- If is unlikely that this is 
correct—a similar “averaging effect” which reduces the vN cross section also affects NN 
collisions. For F(j,ope ~ F", one should consider the true to be a slowly increasing, but 
unknown, function of the density. 

The neutrino signal of SN 1987A offers a unique possibility to test for the actual behavior 
of Fo-. The neutrino signal of a SN depends sensitively on the neutrino opacity at high density, 
which is dominated by Since aA drops precipitously if Fg. > T, we in principle have 
an experimental probe of Fo-(p, T). To this end, three of us (KJR) have performed detailed 
numerical calculations of the cooling of newborn neutron stars with modihed axial-vector 
neutrino opacities @ . The result of this study is that reduced neutrino opacities shorten the 
neutrino pulse dramatically and increase the number of events and the event energies in the 
Kamiokande II and IMB detectors. The best £t to the data is achieved if one ignores opacity 
suppression effects entirely. An overall reduction by more than 50% seems to be excluded. 
This is still compatible with the expected reduction of Ca in a bulk nuclear medium which 
entails an opacity suppression by 30-40% 0, but would be incompatible with any signihcant 
reduction due to either “pairing” of nucleon spins by the NN interaction or broadening of 
Similar conclusions follow from Ref. P] where the opacity was indirectly lowered by 
increasing the number of neutrino species. 

Of course, the SN 1987A observations encompass very few events. Probably the most 
reliable measure of the opacities in the deep interior is the neutrino pulse duration as the 
number of events and their energies are largely determined by the spectral features of the 
neutrino flux and thus by the conditions near the neutron star surface. It is worrying that 
a long signal and thus large opacities are mostly supported by the three late Kamiokande 
events which are separated from the hrst eight by a 7 s break. If the late events had some 
other cause besides the Kelvin-Helmholtz cooling such as late-time accretion, possibly trig¬ 
gered by a reverse-shock , then reduced opacities might become more tolerable. 

Barring this possibility we infer that Fq- cannot grow linearly with density but reaches a 
value of at most a few times T in a SN core. Values for Fg. to the right of the hatched band 
in Fig. ^ seem to be excluded. If correct, this result implies that the OPE value for P^ is an 
overestimate. 

Axion Bounds. Previous axion bounds from SN 1987A had been based on the “naive 
emission rate” (dashed line in Fig. and had employed a perturbative estimate of P^. 
During the first few seconds after collapse Po-,ope/F" = 20-50 throughout the SN core. This 
would take one somewhat to the right of the maximum of the solid line in Fig. If our 
“calibration” of P^ is correct, one is somewhat to the left of the maximum. Either way, Qa 
is near its theoretical maximum so that the “naive” Qa is suppressed by about a factor of 
10. Because Qa is proportional to the square of the axion mass, previous limits on nia are 
relaxed by about a factor of 3. 

The often-quoted limit nia ^ 10“^ eV was based on Cn = 0.5 for V = p and n 
In popular axion models these couplings are much smaller, for example ~ 0 in the so- 
called KSVZ model. Therefore, in a detailed review [Q it was shown that translating the 
SN 1987A bound on the aN coupling into a bound on nia really yields < 3xlO“^eV in 
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the KSVZ and other typical models. Together with the above suppression effect we conclude 
that nia ^ 10“^ eV is not excluded. 

Energy Transfer in vN Collisions. In previous SN studies of neutrino transport, uN 
collisions have been taken to be energy conserving 10.0 because S{uj) oc 5(a;) for scattering 
of Cs from single, “heavy” nucleons. The (and 17^, Ur, and T7r) spectra are formed 
at the radius where ue scattering and ^ become inefficient for transferring 

energy between the neutrinos and the medium. However, the must diffuse through a 
considerable overburden of nuclear material before they escape. In general this material 
will have a temperature T that is lower than the temperature Ty that characterizes the 
spectrum. Thus the average energies are usually found to be a factor of 1.3-1.8 larger than 
those of z7g’s, which react via /3-processes and decouple at lower density and temperature. 

A broad S^ioj) such as in Eq. (5) would allow for energy transfer between the neutrinos 
and the medium. Further, even for vN processes the 5 function is only an approximation 
that neglects recoil corrections. Their inclusion would modify the expected z/^, z7^, Vr, and Vr 
spectra. In turn, energy conservation and a reduced energy flux in these flavors would imply 
increased Ve and V^. energy and number fluxes. This is important in planning for detection 
of SN neutrinos at Superkamiokande or SNO, for understanding the explosion mechanism 
of Type II supernovae 00. and for the conjecture that r-process nucleosynthesis occurs 
in the high-entropy medium that surrounds the protoneutron star 0. 

For neutrinos with a black-body spectrum the average recoil energy transfer in vN colli¬ 
sions with nondegenerate nucleons is {Ae)yN ~ 30Ty(T — Ty)/mN 
energy exchange {Ae)yNN in inelastic uNN scatterings is 


13 , and the average 


/o°° dsiE^e de2£2(^2 — — £ 2 ) 

/o“ de,£^,e-^cn - £2) ' 


( 6 ) 


With Eq. (^, the sum rule for and s(x) = 1 we hnd numerically {Ae)^^^ = C iTo- (T— 

Ty){TyT)~^/‘^ where 1 -|- logio(l + ^^r/T) within 15% for 1 < Ty/T < 2 and 0.1 < 

T„/T < 10. For conditions near the neutrinosphere both the vN and vNN processes can 
therefore transfer several MeV of energy to the medium. Even though the energy exchange 
per interaction may be small, a typical neutrino experiences 10-100 interactions during the 
diffusion process through the protoneutron star surface layers. The net effect may lead to 
substantially lower average energies of the z/^’s and z7^’s. Thus, both vN recoil and the 
inelastic vNN modes of energy transfer are potentially important near the neutrinosphere 
where the spectra are formed. 

The same conclusion is reached by comparing the recoil energy transfer rates between 
neutrinos and a medium of nucleons and between neutrinos and a medium of degenerate 
electrons ||12||, 


p o 1 -1/3 


T 

J- 1/ 


10 MeV 


Y 

e 

or. 


( 7 ) 


Here, the effective weak interaction constant Ce for z/e scattering is 0.33 and 0.26 

(z7^^t-). In the neutron star surface layers p ~ lO^^-lO^^g/cm^, and fy drops to 0.05-0.1 


within about a second after stellar core collapse []^. The energy transfer from z/^ to the 
medium is apparently dominated by direct transfer to the nucleons by an order of magnitude 
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or greater even though the calculation behind Eq. (|^) does not include the full effects of 
electron phase-space blocking. 

Summary. For axial-vector current interactions of neutrinos and axions in a SN core 
multiple-scattering effects must not be ignored. Notably, the spin-fluctuation rate cannot 
be simultaneously very small (Ecr -C T), as had been assumed for the neutrino opacities, and 
very large (Fo- S> T), as had been assumed for the axion emissivities. Our study suggests 
that it is in the middle, i.e. of order T, a range where the impact on neutrino opacities is 
relatively mild while the axion emissivity is near its maximum allowed by multiple-scattering 
effects (but still less than previous calculations had indicated). At the present state of the 
art, the neutrino opacities cannot be calculated from Erst principles within controlled ap¬ 
proximations. Therefore, the SN neutrino signal duration and other signal characteristics 
can be predicted only with the accuracy of a dimensional analysis, implying uncertainties 
for models of the //-driven mechanism of Type II SN explosions |T^. Even near the neu- 


trinosphere nonnegligible effects such as the nonconservation of energy in uN collisions are 
important where both the recoil and the bremsstrahlung mode of uN energy transfer need to 
be included. A quantitative understanding of the impact of energy transfer in uN collisions 
on the expected neutrino fluxes and spectra will require detailed numerical simulations. 

Note Added. After this paper had been submitted, a perturbative calculation of the 
reduction of weak scattering rates by nucleon (iso)spin fluctuation has appeared ||^. In 
essence, this is a calculation of the slope of our (ua) curve in Fig. || at the point T„/T = 0. 
In linear-response theory, this slope is independent of our Lorentzian ansatz and agrees 
with the perturbative result—for a detailed comparison see Ref. |^. Further, Sigl has 
derived an f-sum rule for ^^(a;) which allows one to estimate an upper limit on Tg-, providing 
a theoretical underpinning for our SN 1987A bound on T^/T. 
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FIGURES 



FIG. 1. (a) Bremsstrahlung emission of neutrino pairs in NN collisions. Instead of a neutrino 
pair, an axion can be emitted, (b) “Inelastic vN scattering” in the presence of a bystander nucleon. 



G/r 


FIG. 2. Schematic dependence of the average axial-vector I'N scattering cross section {(ta) and 
of the axion emission rate per nucleon Qalns on T„. The broken lines correspond to the naive 
results which are characterized by F = 0 in Eq. @ for QajuB and by Sa{uj) = 2'k5{uj) for (ua). 












